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ABSTRACT 


New first- and second-order displacement transfer functions have been developed for deformed shape 
calculations of nonuniform cross-sectional beam structures such as aircraft wings. The displacement 
transfer functions are expressed explicitly in terms of beam geometrical parameters and surface strains 
(uniaxial bending strains) obtained at equally spaced strain stations along the surface of the beam 
structure. By inputting the measured or analytically calculated surface strains into the displacement 
transfer functions, one could calculate local slopes, deflections, and cross-sectional twist angles of the 
nonuniform beam structure for mapping the overall structural deformed shapes for visual display. The 
accuracy of deformed shape calculations by the first- and second-order displacement transfer functions 
are determined by comparing these values to the analytically predicted values obtained from finite- 
element analyses. This comparison shows that the new displacement transfer functions could quite 
accurately calculate the deformed shapes of tapered cantilever tubular beams with different tapered 
angles. The accuracy of the present displacement transfer functions also are compared to those of the 
previously developed displacement transfer functions. 

NOMENCLATURE 

c depth factor (vertical distance from neutral axis to bottom surface of uniform beam), in. 
c/ depth factor (vertical distance from neutral axis to z'-th strain station on bottom surface of 
nonuniform beam), in. 

c n value of c; at free end (beam tip), x = x n = / , in. 
c n / Cq beam depth ratio, dimensionless 

c(x) depth factor (vertical distance from neutral axis to bottom surface of nonuniform beam), in. 
cq value of q at fixed end (beam root), x = xq = 0 , in. 
deg degree 

E Young’s modulus, lb/in 2 

i =0, 1, 2, 3, ..., n, strain station identification number 

/ moment of inertia, in 4 

j index 

/ length of beam, in. 

Mi = M(xj ) , bending moment at strain station, X [ , in-lb 
M(x) bending moment at axial location of beam, x , in-lb 
n index for last span-wise strain station (or domain density) 

P force, lb 

SPAR Stuctural Performance And Resizing 
t wall thickness of tubular beam, in. 

x, y Cartesian coordinates (x in axial direction of beam, y in vertical direction), in. 

Xj axial coordinate associated with i -th surface strain station, in. 

yi = y(Xj ) , beam deflection in y-di recti on at axial location, x = Xj , in. 
y(x) beam deflection in y-direction at axial location, x , in. 

a = tan -1 [(cq -c n )l l ] taper angle of beam, deg 

A l = (x; - x,-_i) = l In, distance between two adjacent strain stations { x/_] , x - }, in. 

£j = £(Xj ) , surface bending strain at strain station, x • , in/in 

g(x) surface bending strain (axial strain) at axial location, x , in/in 
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Oj = 9(xj ) , slope of deformed beam at axial location, x = x,- , rad or deg 
0(x) slope of deformed beam at axial location, x , rad or deg 
<7j = E£i , axial stress associated with surface bending strain, £; , lb/in 2 

<t(x) = Ee(x ) , axial stress associated with surface bending strain, e(x), lb/in 2 

INTRODUCTION 

Recently, several displacement transfer functions were formulated for structural deformed shape 
calculations using surface strains (refs. 1-4). In the formulations of the displacement transfer functions, 
the beam-like structure (such as an aircraft wing) was discretized evenly into multiple small domains of 
the same length. Through such discretization, linear-variation assumptions could be made for both beam 
depth and surface strain within each small domain. This approach enabled the integrations of the classical 
beam differential equation in closed forms to yield beam slopes and deflections for each small domain. 
By combining the slope and deflection equations, one can write the final forms of the deflection 
equations, called displacement transfer functions, in terms of beam geometrical parameters and surface 
strains evaluated at evenly distributed inboard strain stations (at domain junctures). Thus, the 
displacement transfer functions can be used to convert the surface strains into displacements (deflections 
and angular rotations) at multiple domain junctures so that one can map out overall deformed shapes of 
the beam structures for visual display. Thus, the displacement transfer functions, with an accompanying 
surface strain-sensing system, has created a revolutionary new structural shape-sensing technology 
(ref. 5, Method for Real-Time Structure Shape-Sensing, Ko-Richards, U.S. Patent No. 7,520,176) for in- 
flight deformed shape monitoring of flexible wings and tails. In addition, the wing shape monitored in 
real time could then be input to the aircraft control system for aeroelastic wing shape control. Similar to 
the structural deformed shape calculations, the distributed surface strains also can be used to predict the 
structural operational loads (ref. 6, Process for Using Surface Strain Measurements to Obtain 
Operational Loads for Complex Structures, Richards-Ko, U.S. Patent No. 7,715,994) for flight vehicle 
operational load monitoring. 

When the structural shape-sensing technology is applied, using the conventional strain gage system 
for sensing the surface strains of flight vehicles is impractical because of excessive lead wire weight. The 
most attractive candidate for a flight vehicle surface strain-sensing system, an alternative to the 
conventional strain gage system, is the fiber-optic strain-sensing system, because the fiber optics are 
lightweight, fine, and flexible filaments (approximately the size of human hairs), and they can be highly 
multiplexed (Bragg gratings) to define the strain-sensing points at desired sensing intervals (refs. 7-9). 
Another powerful characteristic of the fiber-optic strain-sensing system is that, for sensing a given 
structure, the needed number of actual strain data extraction points can be easily specified by means of a 
single command. 

The displacement transfer functions were originally developed for deformed shape predictions of 
straight beams and were successfully validated for accuracy by finite-element analyses of different 
sample structures, such as cantilever tubular beams (uniform, tapered, slightly tapered, stepwise tapered), 
two-point supported tapered tubular beams, flat panels, and tapered wing boxes (unswept and swept). The 
straight-beam displacement transfer functions were further extended to the shape predictions of slender 
curved structures by introducing curvature- effect collections (ref. 10). 

In this report, the first- and second-order displacement transfer functions are developed for deformed 
shape predictions of nonuniform cantilever beam structures (such as aircraft wings). Analytically 
calculated surface strains are then input to the transfer functions, and structural shape calculations are 
made. The shape calculation accuracy of the newly formulated displacement transfer functions are then 
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compared with those of the earlier displacement transfer functions (ref. 2) developed for nonuniform and 
slightly nonuniform cantilever structures. 

BASICS OF DISPLACEMENT TRANSFER FUNCTIONS 

The basics of the displacement transfer functions previously formulated (refs. 1-4) for shape 
predictions of cantilever straight beams using distributed surface strains are briefly described in this 
section. The piecewise-linear assumptions made in the formulation of the displacement transfer functions 
also are discussed. 


Beam Differential Equations 

The formulation of the displacement transfer functions stems from the classical beam differential 
equation (elastic curvature of the deformed straight uniform beam) given by (refs. 1 1-12) 


d 2 y 

dx 2 


M(x) 

El 


( 1 ) 


in which y is the vertical deflection, x is the beam axial coordinate, M(x ) is the bending moment, E is 
the Young’s modulus, and I is the cross-sectional moment of inertia. Equation (1) can be applied to a 
nonuniform beam segment with sufficient accuracy if the beam cross sections change gradually (weak 
nonuniform beam) (ref. 1 1, p. 143). For the weak nonuniform beam case, I is a weak function of x , and 
over a short beam segment, the average value can be used. 

At the cross-section, x, the bending moment, M (x ) , and associated bending strain (axial strain), 
£(x)[= <J(x) / E\, at the bottom surface point of the nonuniform straight beam can be written as 


M(x) = I 


g(x) 

c(x) 


= EI 


£(*) 

c(x) 


( 2 ) 


in which c(x) is the beam depth factor (vertical distance from the neutral axis to the bottom surface of 
the nonuniform beam, or half depth if the neutral axis is located at the half depth of the beam). For a 
tubular beam case, c(x) will be the outer radius. The sign convention of equation (2) is for the positive 
bending moment, M(x) (bending the cantilever beam upward; fig. 1), to induce a positive bending strain, 
£(x), on the bottom surface of the beam. Combining equations (1) and (2), one obtains the modified 
differential equation for the nonuniform beam as 

4 = — ( 3 ) 

dx 2 c(x) 

Note that the modified beam differential equation (3) contains only the beam depth factor, c(x) , and 
the surface bending strain, £(x) . The flexural rigidity, El, is eliminated. The formulation of the 
displacement transfer functions for the nonuniform straight beam is based upon the modified beam 
differential equation (3), which could be sufficiently accurate if the nonuniform beam cross sections 
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change gradually (ref. 11, p. 143). If the mathematical functional forms of (c(x),£(x)} are known, 
equation (3) can be integrated once to yield the beam slope, dy / dx[= tan 9(x )\ , and then integrated a 
second time to yield the beam deflection, y(x) , for each small domain needed for the structural deformed 
shape predictions. 


Discretization 

To carry out the single and double integrations of equation (3) in closed forms, the nonuniform 
cantilever beam of length l was discretized into n number of small domains of equal length, A /(=//«) 
(fig. 1). The z'-th (z' = 0, 1, 2, 3, ..., n) strain stations are to be located on the bottom surface of the 
domain junctures, x = x, (z' = 0, 1, 2, 3, ..., n) (called strain stations, x;). Note from figure 1 that the 
first and last strain stations { Xq , x n } are located at the fixed end (x = xq = 0) and at the free end 
(x = x n = /) , respectively. Within each small domain, x ( -_ ] < x < x ( - , between the two adjacent strain 
stations, { x ; _i ,x - }, the following piecewise linear assumptions can be made with sufficient accuracy for 
slowly changing functions. 


Basic Assumptions 

The following two basic assumptions are needed for deriving the displacement transfer functions: 

1. Within the small domain, x-_j <x <x ; - , between the two adjacent strain stations, { X;_i,x ; - }, the 
depth factor, c(x) , can be represented with a linear function of (x - x,-_] ) , as 

c(x) = c i _ ] -(c i _ ] -c i ) X *j~ ] ; x._, <x<x. (4) 

in which { c;_ \ , q } are the values of c(x) at the strain stations, { X;_j , x- } , respectively. 

2. Within the small domain, x ; -_] < x < x,- , between the two adjacent strain stations, { x,-_j , x,- }, the 
bottom surface bending strain, £(x) , also can be assumed to vary linearly with (x - x ( -_] ) , as 

£(x) = £;•_] - ( e ;_| - £j ) X ; X i~ 1 - x “ x i (5) 

in which { £/_] , £/ } are the values of £(x) at the strain stations, { x ; -_j , x / | , respectively. 

If £(x) is a strong nonlinear function ofx, the domain length, A/, should be reduced so that the piecewise 
linear assumption of equation (5) can remain as a good approximation. 
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Linearity of Bending Moment 

Substituting equations (4) and (5) into equation (2), and expanding the factor, 1 / c(x) , in binomial 
series form because of a small coefficient (q — q_] ) of the slope term (weak nonuniform beam), one 
obtains 


M(x) = El 


£(x) 

c(x) 


= El 


£ i - 1 
c i - 1 


1- 


1--^- 


e i-U 


x-x 


i— 1 


A/ 


1 + 


' c t ' 

1 l - 

v c ;-i y 


X-X/-1 


A/ 


+ 


/ A2 

1 Cj - 

V C i - 1 J 


— — + 

(A/) 2 


( 6 ) 


For a weak nonuniform beam, (1 - q / q_i) is small (that is, (1 - q / q_] ) « 1 ), and within the small 
domain, x;_] < x < x,- , if the average value of / is used, the bending moment, M(x) , within the domain, 
x ; -_i < x < Xj , will vary almost linearly because higher order terms in (1 - q / q_] ) can be neglected. For 
a moderate nonuniform beam, however, the value of (1 — q / q_ j) is no longer very small, and the higher 
order terms must be retained. This retention of the higher order terms will cause M (x) to be slightly 
nonlinear over the domain, x ; -_i < x < x,- . 

For the uniform beam ( c(x) = c = constant), the bending moment, M(x), is directly proportional to 
the bending strain, £/(x) , according to equation (2). Therefore, if the bending moment, M(x) , is a linear 
function of x , then the bending strain, e(x) , also will be a linear function of x . 


Basic Slope Equations 


The slope, tan 6(x)(= dy / dx) , of the nonuniform beam at the axial location, x, within the domain, 
X;_i < x < x; (fig. 1), can be obtained by integrating equation (3) once, and enforcing the continuity of 
the slope at the inboard strain station, x;_] , as 


tan 0(x) = ! 

J X: 


,2 

a y 


dx 


, x i~ 1 dx L 
Integration ofeq. (3) 


+ , an = P ^ 

X-t J i,-1 n(x) 

Slope at x,-_] 


Slope increment 
above tan0 ; _j 


r/x + tan 0,_ ] 
Slope at x,-_i 


( x M < x < x ( - ) 


(7) 


in which tan0/_i is the slope at the inboard strain station, x;_ j. In light of the piecewise linear 
assumptions of (c(x),£(x)} given by equations (4) and (5), respectively, equation (7) can be written in 
the form 


X-Xj - 1 


tan 


£/_i —(£;_] — £;) 

J T_i 


c i — 1 1 c i ) 


A/ 


dx + tan 0 ; _] ; ( x/_ ] < x < x/ ) 


( 8 ) 
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Note that the direct integration of equation (8) will result in a slope equation of the form (ref. 1, eq. 
(23)) 


tan 6 1 = A 1 


£ i - 1 £ i . £ i-\ c i £ i c i - 1 i 
c i-\~ c i (Ci-l-q) 2 


c i 

c i - 1 


+ tan 0/_] 


(1 = 1, 2, 3, n) 


(9) 


in which tan (9/ [ = tan0(x,)] is the slope at the strain station, x ; - . Note from equation (9) that in the limit, 
C'i / c,-_] =1, for the uniform beam, the first term, (£;_] -£/)/(C;_] -c ( -), will go to infinity, and the 

second term with the factor, [log e (c ; - /q_q)]/(q_i -c ; ) , will become mathematically indeterminate 
(that is, 0/0). To circumvent this limit-case mathematical breakdown problem, the logarithmic term must 
be expanded in terms of a (c[_\ — Cj) series in the vicinity of C( / C(_\ ~ 1 , so that the denominators 
containing the factor, (c(_\ - Cj ) , can be cancelled out. After mathematical manipulations are carried out, 
equation (9) can be reduced to the form (ref. 1, eq. (27)) 


n A/ 
tan 0: = 

' 2c. 


i—l 


IY > 

“ 

C: 


2 1 

£ i-l + £ i 

.1 c i- 1 y 



-t-tan#^ ; (/ = 1, 2, 3, ..., n) 


(10) 


which is now applicable to the limit case, q / q_i = 1 , of the uniform beam. 


Basic Deflection Equations 

The deflection, y(x) , of the nonuniform beam within the domain, x;_] < x < x; (fig. 1), can be 
obtained by integrating the slope equation (7) and enforcing the continuity of the deflection at the inboard 
adjacent strain station, x ( -_ i , as 


' — ■ — — ■ — - 1 Deflection ' * ' ' — — — v- — — — ' Deflection 

Integration of slope a tx- Deflection increment above y- , Deflection at x a tx- 

due to tan0-_j 

(x/—] <X<Xj) (11) 

in which y ( -_] is the deflection at the inboard strain station, x/_] . 


rX rX 

y(x)= I tan 0{x)dx + y^_ j = j 

■ x i-\ Jx i - 1 


f 


£(x) 


x i - 1 c(x) 


dx dx 


t-1 


tan G-_ | dx + 


yt - 1 


In light of the piecewise linear assumptions of (c(x),£(x)} given by equations (4) and (5), 
respectively, the deflection equation (11) takes on the form 


1- 

rx rx 

r £ . ^ 

i ' 

i 

1 

l £ i - 1 7 

' A 1 C x 

CIYCI V -1- 1 ton fJ . 1 rj Y V- 1 

*U-1 J x,-_! ^ 

y a 

1 Cf 

liAtiA 1 1 Idll (7^ | liA 1 f j 

X-X,-_l 

v c i-l y 

M 


( 12 ) 
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Note that direct integration of equation (12) will yield the deflection equation in the form (ref. 1, 
eq. (25)) 


>7 = (A/) 


2 £ i - 1 £ i £ i-\ c i £ i c i — 1 


2(q_i-q) (c M - C/ ) 3 


Ci log £ 


c; 


C; 


- + (q_ | - q ) 


/-I 


• + y i _l+A/tan0 ; _ 1 

(7 = 1, 2, 3,..., n) (13) 


in which y* [= y(x; )] is the deflection at the strain station, x- . Note from equation (13) that in the limit, 
C{ / C(_\ = 1 , for a uniform beam case, the first and third terms with the factor, (q_q-q), in the 

denominators will go to infinity, and the second term with the factor, [log e (c ; - / Cf_\ )]/('e ( -_| — C() , will 
become mathematically indeterminate (that is, 0/0). To circumvent this limit-case mathematical 
breakdown problem, the logarithmic term must be expanded in terms of a (c ; -_] — C[) series in the 
neighborhood of c; / c ( -_] ~ 1 so that the denominators containing the factor, (c ; -_] — c/ ) , can be cancelled 
out. After the mathematical manipulations are carried out, equation (13) can be reduced to the form 
(ref. 1, eq. (30)) 


yi 


(A If 


6c, 


7—1 


IY A 


3- C ' 

e 7-l + e 7 

_V c i - 1 j 



+ yi-\ + A/tan0 ; _i ; (7 = 1, 2, 3, ..., n) 


(14) 


which is now applicable to the limit case, c; / = 1 , of the uniform beam. 

ALTERNATIVE DISPLACEMENT TRANSFER FUNCTIONS 

To circumvent the aforementioned limit-case mathematical problems, the development of alternative 
displacement transfer functions for the nonuniform cantilever beam is needed. The new alternative 
displacement transfer functions also must be applicable to the limit case of a uniform beam without 
creating limit-case mathematical indeterminacy problems. Because the cantilever beam is assumed to be a 
weak nonuniform beam (that is, 1 — (c; / c ; -_] ) « 1 ), the perturbation method can be applied to formulate 
the new displacement transfer functions by expanding the factor, 1 / c(x) , in the integrand of equations (8) 
or (12) in binomial series form up to the first- or second-order terms. The formulations of the new first- 
and second-order displacement transfer functions are presented in the following sections. Appendices A 
and B present the detailed mathematical derivations for readers to easily follow. 

FIRST-ORDER DISPLACEMENT TRANSFER FUNCTION 

In this section, the first-order slope equations and first-order deflection equations are formulated for 
the nonuniform cantilever beam. Appendix A presents the detailed mathematical derivations. 
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First-Order Slope Equations 


Under the assumption of the present weak nonuniform cantilever beam for which (1 - c ; - / C (_\ ) « 1 , 


the factor, 


i-l 


1- 


l--^- 
( c i- 1 ) 


X-Xi-l 


A / 


in the integrand of the slope equation (8) can be expanded in 


binomial series form up to the first-order term to yield 


tan 


0 (*) = ^P 

c. , J x- 


i-l "/-I 


1- 


1--^- 


V e i- 1 ) 


x — x 


i-l 


M 


1 + 


1--^- 


x-x 


i-l 


v L i-i ; 


A/ 


+ ... 


dx + tan 0-_ | 


(x-_!<x<x.) (15) 


in which tan 0-_ | f= tan 0(x/_] )1 is the slope at x = x ( -_| . 


After carrying out the integration of equation (15) (ref. 13), one obtains the first-order slope equation 
for the domain, x ;_ \ < x < Xj (see Appendix A for detailed derivations), 


tan 0(x) = 


e i-l 


c; 


i— 1 


(X — X;_l) + 


C; 


V £ i— 1 c i— 1 


(x-Xi_\) z 


2M 


1 

V £ i-\ 


c; 


1- 

v c i-i ; 


( x ~ x i-V 


3(A lY 


+ tan Q;_ 


'i-l 


(x-_| <x<x- ) (16) 


At the strain station, X(, we have (x ; - - X;_i ) = A/ , and equation (16) takes on the form (see 
Appendix A) 


tan 6/ = 


A 1 


6c, 


i — 1 


r r . \ ( 

4 l — |£;_| + 


V 


C; 


i-l y 


5-2- 


V 


Ci 


i-l ) 


+ tan Oi_i ; (i = l, 2, 3, ..., n) 


(17) 


in which tan 0/ [ = tan d(xj )1 is the slope at x = x ; - . Equation (17) is the first-order slope equation in a 
descending recursion relationship. Applying the descending recursion relationships, one can rewrite 
equation (17) in a series summation form for the strain station, x ( - , as 


^ A/ v 1 

tan0 i = "T X 

6 . c 



y 



7 

X 

“ 


4 — 

C J 

£ /'-l + 

5- 

V 

c i 
0 J 

£ ./ 

■1 

A 

c 7-lJ 

Z. 

c 7-lJ 


+ tan0p ; 0 = 1, 2, 3, 

=0 for canti- 
lever beams 


«) 


( 18 ) 



Equation (18) is the final form of the first-order slope equation in series summation form for the 
nonuniform cantilever beam. 

For the uniform beam case ( c ■ = c-_ j = c ), equation (18) degenerates into the form 


A / 


tandi = — ^(£j-i+£j)+ tan0 o 


2c 


0 = 1, 2, 3, n) 


7=1 


=0 for canti- 
lever beams 


(19) 


which agrees with previously established equation (25) of reference 2. 


First-Order Deflection Equations 

The deflection, y(x ) , of the nonuniform cantilever beam within the domain, jc /_ \ <x<x ; -, can be 
obtained by integrating the slope equation (16) (ref. 13), and enforcing the continuity of the deflection at 
the inboard adjacent strain station, x ( -_] , as (see Appendix A) 
rx 


pX 

y(x)= tan 9(x)dx + >’/_] 
' “ eq. (16) 


(x-Xj_if £/- 1 


12 


6 + 2 


c i-l 


c; 


V fy — 1 c /-l y 


(x-x,_i) 


A/ 


1--SL- 


£f- 1 


1- 


C; 


c i-l y 1 


(X-X;_|)" 


(A/r 


+ (x-x f _ 1 )tan0 f _ 1 +y i _ 1 

( ^/- l < x < x,- ) 


( 20 ) 


in which y ; -_] is the deflection at the inboard strain station, x/_] . 


At the strain station, x, , we have (x,- - x/_] ) = A/ , and equation (20) yields the deflection in the form 
(see Appendix A) 


y'i 


(AlY 

12c 


7—1 




V 


c; 


i-l 7 


3 — — 


V 


C; 


l-l 7 


+ yi-\ + A/tan0 ; _i ; (i = l, 2, 3, ..., «) 


( 21 ) 


in which yi[=y(xj )] is the deflection at the strain station, x ; - . Equation (21) is the first-order deflection 
equation in a descending recursion relationship. 
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For a uniform beam (q = q_] = c), equation (21) degenerates into the form 


yt = 


(A/) 2 

6c 


(2e,-_i + £/) + >V-l + A/ tan 0/_] 


; (7 = 1, 2, 3, n) 


( 22 ) 


which agrees with the previously established equation (26) of reference 2. 

Substituting the slope equation (18) into the deflection equation (21), with descending recursion 
relationships applied, yields (see Appendix A for derivations) 


2 7 


yt 


(A/) 


12 (L 


x- 


5 + 8(;-l)-(l + 2(;-l)>^±i 


c i~J 


£ i-.i 


+ 


3 + 1 0(7 - 1) - (1 + 4(7 - 1)) - 


c i~j + 1 
c i~j 


H-j + 1 f + yp + (i)Al tan flp 

=0 for cantilever beams 
(i=l, 2, 3, ..., n) 


(23) 


Equation (23) is the final form of the first-order deflection equation in series summation form (called the 
first-order displacement transfer function) for the nonuniform cantilever beam. 

For the uniform beam case ( q = q_i = c ), equation (23) degenerates into 

( A l 

yi = X [cv - l ) £ i-j + Oj - 2 ) £ i-j +\ ] + )’() + (i)AI tan 6 () ; (i = 1, 2, 3, ..., n) (24) 

J = 1 =0 for cantilever beams 

which is identical to equation (27) of reference 2. 

After the terms are grouped, equation (24) can be rewritten in the alternative form (see Appendix A 
for derivations) 


>7 = 


(A/) 2 
6 c 


i — 1 

(3/ - 1)£q + 6 ^ (i - j)£ j + £[ 
7=1 


+ )'o + (i)AI tan d() 

=0 for cantilever beams 


(7 = 1, 2, 3, ..., n ) 


(25) 


which agrees with equation (28) of reference 2. 
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SECOND-ORDER DISPLACEMENT TRANSFER FUNCTION 


In this section, the second-order slope equations and second-order deflection equations are formulated 
for the nonuniform cantilever beam. Appendix B presents the detailed mathematical derivations. 

Second-Order Slope Equations 

For the weak nonuniform beam with a slowly changing depth (that is, (1 - c ; - / c ; -_] ) « 1 ), accuracy of 


the displacement transfer functions can be improved by expanding the factor, 


i—l 


1- 


l--^- 
V c i- 1 ) 


x-Xi-l 


A / 


in the integrand of the slope equation (8) in binomial series form up to the second-order term as 


tan 0(x) = 


£ i - 1 ( x 


J 

J X: 


Ci-\ JX i - 1 


1- 


1--^- 


£ i - 1 


X-Xi_ 1 


A / 


1 + 


l--^- 
Ci - 1 




A/ 


1--3- 

Q-l 


(x - x ; -_| 


+ .... 


rfx 


+ tan0,_ 


7-1 


(A IT 

(x._i<x<x.) (26) 


in which tan0,_i is the slope at the inboard strain station, x/_] . 

After carrying out the integration of equation (26), one obtains the second-order slope equation for 
the domain, x/_ ] < x < X [ , as (see Appendix B for details) 


tan 0(x) = 


gf-1 
c i - 1 


(x-x,-_i) + 


c,- 


V £ i — 1 c i — 1 


(x-x,_i) 


2A/ 


t c i J7 (7_ 

l C i— 1 y V — 1 C i— 1 


(X-X;_l) J 


3(A iy 


( \ 2 f 

V c i- 1 ; 


1--* 


l gf-l 


(X-X;_|) 


4(A/) J 


+ tan 0, 


7-1 


(x,_i <x<x f ) (27) 


At the strain station, x = x ( - , we have (x/ - x/_] ) = A/ , and equation (27) becomes (see Appendix B) 


tan 0/ 


V 


12c,-_ 


7—1 




9 — 4- 


<7-1 


i c t-i ; 


£,•-1 + 


13-10 — — + 3 
c i - 1 


6 3 

2' 


c i 





gf 

V <7-1 y 




+ tan0 ; _i 

(i' = l, 2, 3, ..., n) (28) 
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in which tan0/[s tan d(x[ )1 is the slope at the strain station, X { . Equation (28) is the second-order slope 
equation in a descending recursion relationship. Applying the descending recursion relationships, one can 
write equation (28) in series summation form as 


n A/ v 1 

tan 6j = — V < 

12 ^ c • i 1 
j = 1 


c i 

9-4 — — + 


f c - 3 

C J 


C J-1 


V C J~ l J 


£ j~ 1 + 


13-10-^— + 3 


C J-1 


( \ 

2' 


c j 





£ j 

Uj-iJ 




+ tan 0() 

=0 for canti- 
lever beams 

(i=\, 2, 3, ..., n) (29) 


Equation (29) is the final form of the second-order slope equation in series summation form for the 
nonuniform cantilever beam. 

For the uniform beam case (c ; - = Cj_\ = c), the second-order slope equation (29) degenerates into 


Al 1 

tan0/ = 2c^ ( £ j~ l+£ j^ + ^ = ’’ 2 ’ 3 ’ 

j— 1 =0 for canti- 

lever beams 


(30) 


which agrees with equation (25) of reference 2. 

Second-Order Deflection Equations 

The deflection, y(x ) , of the nonuniform beam within the domain, x;_i < x < x,- , can be obtained by 
integrating the slope equation (27), and enforcing the continuity of the deflection at the inboard adjacent 
strain station, x ( _ ] , as 


fX 

y(x) = I tan 6(x)dx + y ; -_] 

1 eq. (27) 

= 7^ (x-X;_i) 2 +^ 

2c i - 1 3 


c; 




V £ i— 1 c i— 1 


(x — Xj_ i ) | 1 
Al 6 




1 -- 


c; 


c; 


1 

To 




\ c i - 1 J 


1--* 


V £ /-i 


(x — Xl_\f 


(Air 


V c i— 1 J\ £ i — 1 c i — 1 

f +»-l +(X-X;_] )tan0;_| 


(x-x,-_i) 


(A/)" 


(X/-1 <x<Xi) (31) 


in which is the deflection at the inboard strain station, x ; -_j . 
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At the strain station, X ( , we have (xj - x,-_] ) = A/ , and equation (31) becomes (see Appendix B) 


Tr 


(a iy 

60 c, 


i - 1 


27-9— 1 — + 2 
c 7-i 


/ A2 

_fi_ 


£ i - 1 + 


18-11— 1 — + 3 
c 7-l 


( \ 

2~ 

' 

Cj 



l 


£j 

y C 7— 1 J 




+ ^._ 1 +A/tan6>._ 1 
(7 = 1, 2, 3, n) 


( 32 ) 


which is the second-order deflection equation in a descending recursion relationship. 
For the uniform beam case (c/ = c,-_] = c), equation (32) degenerates into 


yi 


(Air 

6c 


■(2e/_i + £,-) + y;_i + A/tan0 ; -_i ; (7 = 1, 2, 3, n) 


(33) 


which is identical to equation (26) of reference 2. 

Substituting the slope equation (29) into the deflection equation (32), with descending recursion 
relationships applied, yields (see Appendix B) 


>7 = 


(A/) 


,2 i i 
60 c i~j 



( \ 

2" 

27 + 45(j - 1) - (9 + 20(7 - 1)) + (2 + 5(j - 1)) 

c 7— 7+1 


c i~j 

v Ci ~j > 



£ i~j 


18 + 65(j - 1) - (1 1 + 50(j - 1)) — — ■ + (3 + 1 5(_/ - 1)) 


C i-J 


( \ 

2" 


c 7— 7+1 





£ 7— 7+1 

l c i-j j 




+ )'o + (/)A/ tan 0() 

v ' v ' 

=0 for cantilever beams 

(7 = 1, 2, 3, n) (34) 


Equation (34) is the final form of the second-order deflection equation in series summation form (called 
the second-order displacement transfer function) for the nonuniform cantilever beam. 

For the uniform beam case (c,- = q _ \ = c), equation (34) degenerates into 


>7 = 


(Air 


6c 


X @j - \)£i-j + (3 j - 2 )£i_j+i ] + yp + (i)AI tan 6 () ; (7 = 1, 2, 3, ..., n) 


(35) 


7=1 


=0 for cantilever beams 


which is identical to equation (27) of reference 2. 


SUMMARY OF ALL DISPLACEMENT TRANSFER FUNCTIONS 

This section summarizes all the displacement transfer functions previously developed (refs. 1, 2, 4). 
These transfer functions are listed for comparison of mathematical functional forms. 
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Previous Displacement Transfer Functions 

The displacement transfer functions previously developed for nonuniform and slightly nonuniform 
straight cantilever beams have the following mathematical forms (refs. 1, 2, 4): 

1. The displacement transfer function for the nonuniform cantilever straight beam ( q ^ q_] ) 
(derivations in ref. 2, Appendix A) is expressed as 


y i = (M) 2 J j < 
7=1 


[2(i — j) + 1] 


£ j ~ 1 fj £ /-l c J £ 7 c 7-1 
2(c_/-i-c 7 ) (c m -c 7 ) 3 


c.log 


c y-i 


+ (c,_ 1 -c / ) 


i — 1 

+(A/) 2 £(/-;) 
7=1 


£ 7 ~ 1 c 7 ^ 9 -! 


log 



+ V() + (i)A/ tan 0() 

v v J 

=0 for cantilever beams 


(i = 1, 2, 3, ..., n) (36) 

Equation (36) is not applicable to the uniform beam because of mathematical indeterminacy 
problems at q = c ( -_| = c . 


2. The displacement transfer function for the slightly nonuniform cantilever straight beam 
( C'i / c/_] — > 1 ) is expressed as 


}’i = 


(Air 


j=\ Ci ~J IL 


3(2; — 1) — (3./ — 2)- 


c i~j + 1 
c i~j 


£i-j + (3 j - 2)£j_j +] U y 0 + (i)Al tan 9 0 


=0 for cantilever beams 

(i=\, 2, 3, ..., n) (37) 


which was obtained from equation (36) by expanding the logarithmic terms in series form in the 
neighborhood of q / ~ 1 (derivations in ref. 2, Appendix C). 


New Displacement Transfer Functions 

The two new displacement transfer functions (23) and (34) developed in this report are duplicated as 
equations (38) and (39), respectively, as follows: 

1. The first-order displacement transfer function (eq. (23)) is expressed as 


yi 


(A/) 


2 i 


X- 


12 (L 


5 + 8(; - 1) - <1 + 2(; - 1)) 


c i~j + 1 
c i~j 


£ i~j 


+ 


3 + 10(; - 1) - (1 + 4(; - 1)) - 


c i~j + 1 
c i~j 


H-j+ lf+ TQ + (QA/ tan 6 >q 

=0 for cantilever beams 

(i=\, 2, 3, ..., n) (38) 
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2. The second-order displacement transfer function (eq. (34)) is expressed as 


O i 


( \ 

2’ 

(A/) ^ 1 

yr ~ « h*. / 

7=1 l ~J 

27 + 45(y - 1) - (9 + 20C/ - 1)) + (2 + 5(y - 1)) 

c i-j+l 


C i-j 

\ l ~J J 



' i-J 


18 + 65(7 — 1) — (11 + 50(7 — l))^ ±i + (3 + 15(7 — 1)) 

( \ 

*7—7+1 

2 ’ 

*7 7+1 

c i~j 

V l ~J ) 




+ 


y () + (i)Al tan0Q 
=0 for cantilever beams 

0 = 1, 2, 3, n) (39) 


3. The displacement transfer function for the uniform straight cantilever beam (c;=c/_i = c) 
(eq. (25)) is expressed as 


y,= 


(A/) 2 

6c 


i — 1 

(3/ - 1)£q + 6 ^ (/ - j)£j + £■ 
7=1 


+ jVq + (/)A/tan0Q 
=0 for cantilever beams 


0=1, 2, 3, n ) 


(40) 


Equation (40) is the degenerated form of equations (37), (38), and (39) for the uniform beam case 
(q = q_ i = c) (see Appendices A and B). 


Characteristics of Displacement Transfer Functions 

Except for equation (36), equations (37), (38), and (39) could degenerate into equation (40) for the 
uniform beam case ( c/ = q_] = c) without a mathematical indeterminacy problem (see Appendix A and 

Appendix B). Note that in each of the displacement transfer functions (36)-(40), the deflection, y - , at the 
current strain station, X[ , is expressed in terms of domain length, A/, inboard beam half depths 
( Cq , Cj , c 2 , Cj , . . . , Cj ), and the associated inboard strains ( £q , £^ , £ 0 , £3 ,...,£.), including the values of q 
and £ ; - at the current strain station, X[ , where the deflection, , is calculated. Because deflection 

equations (36)-(40) contain no structural properties, in the deformed shape predictions of complex 
structures (such as aircraft wings), one can avoid tedious computations of bending stiffness, El , at 
different strain stations (cross sections). In fact, the effect of El is absorbed implicitly by the surface 
bending strains, £/ . Thus, not needing to know El is the powerful characteristic of the displacement 
transfer functions for structural deformed shape predictions. 

ANALYTICAL SHAPE PREDICTIONS 

The shape calculations in this report are referred to as an analytical shape prediction study. Namely, 
instead of using actual measured surface bending strains (which require actual flights or time-consuming 
ground tests), the Structural Performance And Resizing (SPAR) finite-element computer program 
(ref. 14) was used to generate the surface bending strains and beam deflections. The SPAR-generated 
deflection curves (assumed to be the correct deflection curves) were then used as reference points to 
check the accuracy of the corresponding deflection curves calculated from different displacement transfer 
functions. 
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The surface bending strains, £ ; - (7 = 0, 1, 2, 3, n ) , at the z'-th strain stations needed for inputs to 
calculate the beam deflections were generated by converting the SPAR-generated axial nodal stresses, 
0 [ (z = 0, 1, 2, 3, ..., n ) , at the z'-th strain station into the associated bending strains, £ ; - , through Hooke’s 

law (that is, £,• = <7/ IE), because the lateral stresses are practically zero for beam structures. Alternatively, 
the surface bending strains, £/ , also can be generated from the axial length changes of the SPAR 
elements where the bending strains, £,• , are to be obtained. This method was found to lose accuracy, 
however, near the highly bent beam tip region and therefore was not used. For accuracy verification, the 
SPAR-generated bending strains were then input to different displacement transfer functions for 
calculations of theoretical deflection curves for comparison with the corresponding SPAR-generated 
deflection curves. 

COMPARISON OF SHAPE PREDICTION ACCURACY 

In the following section, the shape prediction accuracy of the first- and second-order displacement 
transfer functions are compared with those of previous displacement functions developed for nonuniform 
and slightly nonuniform cantilever beams. For the comparative shape prediction accuracy studies, tapered 
cantilever tubular beams with different taper angles were chosen. 

Tapered Cantilever Tubular Beams 

For the tapered cantilever tubular beam (fig.l), the equally spaced strain stations will be along the 
bottom generatrix of the beam. Each of the tapered cantilever tubular beams analyzed has a radius that 
linearly decreases from the fixed end toward the free end. Table 1 lists the dimensions and taper angles, 
a , of the aluminum tapered tubular beams considered. 


Table 1. Dimensions of aluminum tapered cantilever tubular beams. 


/, in. 
(length) 

t, in. 

(wall thickness) 

c Q , in. 
(root depth) 

c n ’ in - 
(tip depth) 

c n ,c 0 

(depth ratio) 

a {=tan _ 1 [(c 0 -c w )//]}, 
deg (taper angle) 

100.5 

0.02296 

4 

4 (uniform) 

4/4 (1.0) 

0.00 

100.5 

0.02296 

4 

3 

3/4 (0.75) 

0.57 

100.5 

0.02296 

4 

2 

2/4 (0.5) 

1.14 

100.5 

0.02296 

4 

1 

1/4 (0.25) 

1.71 

100.5 

0.02296 

4 

0.5 

0.5/4 (0.125) 

1.99 

100.5 

0.02296 

4 

0.25 

0.25/4 (0.0625) 

2.14 


Figure 2 shows the finite-element model of a typical tapered cantilever tubular beam (depth ratio of 
c n /co= 2/4) generated from the SPAR finite-element computer program (ref. 14). The size of the SPAR 
model (identical for all tapered beam cases) also is indicated in figure 2. The tapered cantilever tubular 
beam is subjected to an upward load of P = 100 lb at the beam tip. Figure 3 shows the SPAR-generated 
deformed shapes of the tapered cantilever tubular beams with different taper angles. For the higher taper 
angles (figs. 3d-f), the tapering down region near each beam tip is highly bent, and therefore, high 
bending strains are induced in this region (see next section). 

Surface Bending Strains 

Figure 4 shows the SPAR-generated surface bending strain curves for different tapered cantilever 
tubular beams under bending, using the domain density, n = 8 (the total number of strain stations is 
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\ + n = 9 ). The bending strain curve for the uniform beam case ( c n / c o = 4/4) is a classical straight line. 
For the tapered beams, however, the strain curves are bow shaped, and the bow bend becomes deeper as 
the depth ratio, c h /cq, decreases. At the depth ratio, c n / Cq = 0.5/4, the strain curve is deeply bent by 

nearly 130° toward the beam tip. When the depth ratio reaches the minimum of c n /cq = 0.25/4, the strain 

curve changes to a “cat ear” shape with a nearly 160° sharp bend near the beam tip region where the 
strain reaches a maximum, and then sharply drops toward the beam tip. Note from figure 4 that the 
SPAR-generated strains are not exactly zero (theoretically zero) at the free beam tip (especially for the 
c„/c 0 = 0.5/4 and c n /cQ= 0.25/4 cases) because of the discretization of the beam with finite elements 

(fig- 2). 

Keep in mind that in the formulation of the present first- and second-order displacement transfer 
functions and all the previous displacement transfer functions, the distribution of the bending strain along 
the discretized cantilever beam was assumed to be piecewise linear (eq. (5)). Therefore, the strain curves 
used in the formulation should be piecewise-linear curves (connecting the two adjacent strain data points 
with straight lines) instead of the smooth strain curves shown in figure 4. For moderately bent strain 
curves (fig. 4), the piecewise-linear approximation could be quite accurate, but in a strongly bent region 
of the strain curve (that is, c w /cq= 0.25/4, fig. 4), the piecewise linear approximation may lose some 

accuracy locally using the domain density of n = 8. Thus, for sharply bent strain curves, more strain data 
points (that is, increasing domain density, n) are needed to make the piecewise-linear strain curve 
approach the true smooth strain curves for improving the shape prediction accuracy (see the section 
entitled, “Domain Density”). 


Deflection Data 

The strain data shown in figure 4 were input to deflection equations (36)-(40) for calculations of 
deflections of the tapered cantilever tubular beams with six depth ratios: c n / cq = 4/4, 3/4, 2/4, 1/4, 0.5/4, 

0.25/4. Tables 2a-f compare the resulting predicted deflections with the corresponding SPAR-generated 
deflections. The prediction errors produced by deflection equations (36)-(39) are indicated in parentheses 
for comparison. The percent prediction error at a strain station, jt; , produced by each deflection equation 
is defined as the difference between the predicted and SPAR-generated deflections divided by the 
associated SPAR-generated beam tip deflection. 


Table 2a. Comparison of deflections calculated from SPAR with those calculated from different 
deflection equations; uniform cantilever tubular beam; c n / Cq = 4/4; n = 8. 


Theories 

Jo 

(error, 

percent) 

y i 

(error, 

percent) 

J2 

(error, 

percent) 

J3 

(error, 

percent) 

J4 

(error, 

percent) 

J5 

(error, 

percent) 

J6 

(error, 

percent) 

J7 

(error, 

percent) 

J8 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01787 

(0.0000) 

0.06340 

(0.0000) 

0.13368 

(0.0000) 

0.22450 

(0.0000) 

0.33181 

(0.0000) 

0.45174 

(0.0000) 

0.57933 

(0.0000) 

0.71162 

(0.0000) 

Nonuniform, 
eq. (36) 



















Slightly 
nonuniform, 
eq. (37) 

0.00000 

(0.0000) 

0.01573 

(0.3007) 

0.06027 

(0.4398) 

0.12952 

(0.5846) 

0.21936 

(0.7223) 

0.32568 

(0.8614) 

0.44435 

(1.0385) 

0.57127 

(1.1326) 

0.70230 

(1.3097) 

First-order, 
eq. (38) 

0.00000 

(0.0000) 

0.01573 

(0.3007) 

0.06027 

(0.4398) 

0.12952 

(0.5846) 

0.21936 

(0.7223) 

0.32568 

(0.8614) 

0.44435 

(1.0385) 

0.57127 

(1.1326) 

0.70230 

(1.3097) 

Second-order, 
eq. (39) 

0.00000 

(0.0000) 

0.01573 

(0.3007) 

0.06027 

(0.4398) 

0.12952 

(0.5846) 

0.21936 

(0.7223) 

0.32568 

(0.8614) 

0.44435 

(1.0385) 

0.57127 

(1.1326) 

0.70230 

(1.3097) 


17 




Table 2b. Comparison of deflections calculated from SPAR with those calculated from different 


deflection equations; tapered cantilever tubular beam; c n /cq = 3/4; n = 8. 


Theories 

Jo 

(error, 

percent) 

Jl 

(error, 

percent) 

J2 

(error, 

percent) 

J3 

(error, 

percent) 

J 4 

(error, 

percent) 

y 5 

(error, 

percent) 

J6 

(error, 

percent) 

J? 

(error, 

percent) 

Jg 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01823 

(0.0000) 

0.06697 

(0.0000) 

0.14572 

(0.0000) 

0.25204 

(0.0000) 

0.38294 

(0.0000) 

0.53428 

(0.0000) 

0.70092 

(0.0000) 

0.87634 

(0.0000) 

Nonuniform, 
eq. (36) 

0.00000 

(0.0000) 

0.01623 

(0.2282) 

0.06407 

(0.3309) 

0.14180 

(0.4473) 

0.24707 

(0.5671) 

0.37676 

(0.7052) 

0.52686 

(0.8467) 

0.69210 

(1.0065) 

0.86568 

(1.2164) 

Slightly 
nonuniform, 
eq. (37) 

0.00000 

(0.0000) 

0.01624 

(0.2271) 

0.06409 

(0.3286) 

0.14185 

(0.4416) 

0.24718 

(0.5546) 

0.37699 

(0.6790) 

0.52724 

(0.8033) 

0.69269 

(0.9391) 

0.86659 

(1.1126) 

First-order, 
eq. (38) 

0.00000 

(0.0000) 

0.01623 

(0.2282) 

0.06405 

(0.3332) 

0.14175 

(0.4530) 

0.24699 

(0.5763) 

0.37665 

(0.7178) 

0.52669 

(0.8661) 

0.69187 

(1.0327) 

0.86538 

(1.2507) 

Second- 
order, 
eq. (39) 

0.00000 

(0.0000) 

0.01623 

(0.2282) 

0.06407 

(0.3309) 

0.14179 

(0.4485) 

0.24706 

(0.5683) 

0.37676 

(0.7052) 

0.52686 

(0.8467) 

0.69209 

(1.0076) 

0.86567 

(1.2176) 


Table 2c. Comparison of deflections calculated from SPAR with those calculated from different 


deflection equations; tapered cantilever tubular beam; c n / cq = 2/4; n = 8. 


Theories 

Jo 

(error, 

percent) 

Jl 

(error, 

percent) 

J2 

(error, 

percent) 

J3 

(error, 

percent) 

J4 

(error, 

percent) 

J5 

(error, 

percent) 

J6 

(error, 

percent) 

J7 

(error, 

percent) 

Jg 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01864 

(0.0000) 

0.07113 

(0.0000) 

0.16057 

(0.0000) 

0.28831 

(0.0000) 

0.45510 

(0.0000) 

0.65965 

(0.0000) 

0.89759 

(0.0000) 

1.15801 

(0.0000) 

Nonuniform, 
eq. (36) 

0.00000 

(0.0000) 

0.01681 

(0.1580) 

0.06855 

(0.2228) 

0.15709 

(0.3005) 

0.28385 

(0.3851) 

0.44936 

(0.4957) 

0.65236 

(0.6295) 

0.88815 

(0.8152) 

1.14592 

(1.0440) 

Slightly 
nonuniform, 
eq. (37) 

0.00000 

(0.0000) 

0.01680 

(0.1589) 

0.06849 

(0.2280) 

0.15693 

(0.3143) 

0.28357 

(0.4093) 

0.44895 

(0.5311) 

0.65188 

(0.6710) 

0.88778 

(0.8471) 

1.14549 

(1.0812) 

First-order, 
eq. (38) 

0.00000 

(0.0000) 

0.01679 

(0.1598) 

0.06848 

(0.2288) 

0.15690 

(0.3169) 

0.28348 

(0.4171) 

0.44873 

(0.5510) 

0.65137 

(0.7150) 

0.88672 

(0.9387) 

1.14334 

(1.2668) 

Second- 
order, 
eq. (39) 

0.00000 

(0.0000) 

0.01681 

(0.1580) 

0.06855 

(0.2228) 

0.15708 

(0.3014) 

0.28383 

(0.3869) 

0.44932 

(0.4991) 

0.65230 

(0.6347) 

0.88807 

(0.8221) 

1.14518 

(1.1079) 


Table 2d. Comparison of deflections calculated from SPAR with those calculated from different 


deflection equations; tapered cantilever tubular beam; c n / cq = 1/4; n = 8. 


Theories 

Jo 

(error, 

percent) 

Jl 

(error, 

percent) 

J2 

(error, 

percent) 

J3 

(error, 

percent) 

J4 

(error, 

percent) 

J5 

(error, 

percent) 

J6 

(error, 

percent) 

J7 

(error, 

percent) 

Jg 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01912 

(0.0000) 

0.07598 

(0.0000) 

0.17923 

(0.0000) 

0.33816 

(0.0000) 

0.56506 

(0.0000) 

0.87454 

(0.0000) 

1.28196 

(0.0000) 

1.78417 

(0.0000) 

Nonuniform, 
eq. (36) 

0.00000 

(0.0000) 

0.01742 

(0.0953) 

0.07372 

(0.1267) 

0.17627 

(0.1659) 

0.33448 

(0.2063) 

0.56030 

(0.2668) 

0.86823 

(0.3537) 

1.27201 

(0.5577) 

1.76228 

(1.2269) 

Slightly 
nonuniform, 
Eq. (37) 

0.00000 

(0.0000) 

0.01738 

(0.0975) 

0.07344 

(0.1424) 

0.17549 

(0.2096) 

0.33280 

(0.3004) 

0.55718 

(0.4417) 

0.86290 

(0.6524) 

1.26376 

(1.0201) 

1.75295 

(1.7498) 

First-order, 
eq. (38) 

0.00000 

(0.0000) 

0.01739 

(0.0970) 

0.07355 

(0.1362) 

0.17577 

(0.1939) 

0.33339 

(0.2674) 

0.55820 

(0.3845) 

0.86444 

(0.5661) 

1.26542 

(0.9270) 

1.75143 

(1.8350) 

Second- 
order, 
eq. (39) 

0.00000 

(0.0000) 

0.01742 

(0.0953) 

0.07371 

(0.1272) 

0.17623 

(0.1681) 

0.33440 

(0.2107) 

0.56013 

(0.2763) 

0.86788 

(0.3733) 

1.27135 

(0.5947) 

1.76105 

(1.2958) 
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Table 2e. Comparison of deflections calculated from SPAR with those calculated from different 


deflection equations; tapered cantilever tubular beam; c n /cq = 0.5/4; n = 8. 


Theories 

Jo 

(error, 

percent) 

Jl 

(error, 

percent) 

J2 

(error, 

percent) 

J3 

(error, 

percent) 

y a 

(error, 

percent) 

Js 

(error, 

percent) 

J6 

(error, 

percent) 

J? 

(error, 

percent) 

Jg 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01938 

(0.0000) 

0.07872 

(0.0000) 

0.19047 

(0.0000) 

0.37082 

(0.0000) 

0.64535 

(0.0000) 

1.05517 

(0.0000) 

1.67290 

(0.0000) 

2.59712 

(0.0000) 

Nonuniform, 
eq. (36) 

0.00000 

(0.0000) 

0.01775 

(0.0627) 

0.07665 

(0.0797) 

0.18790 

(0.0990) 

0.36790 

(0.1124) 

0.64197 

(0.1371) 

1.05147 

(0.1425) 

1.66633 

(0.2530) 

2.54871 

(1.8640) 

Slightly 
nonuniform, 
eq. (37) 

0.00000 

(0.0000) 

0.01768 

(0.0656) 

0.07620 

(0.0970) 

0.18656 

(0.1506) 

0.36482 

(0.2310) 

0.63559 

(0.3758) 

1.03867 

(0.6353) 

1.64049 

(1.2479) 

2.50844 

(3.4146) 

First-order, 
eq. (38) 

0.00000 

(0.0000) 

0.01771 

(0.0643) 

0.07640 

(0.0893) 

0.18714 

(0.1282) 

0.36615 

(0.1798) 

0.63829 

(0.2718) 

1.04386 

(0.4355) 

1.64986 

(0.8871) 

2.51231 

(3.2655) 

Second- 
order, 
eq. (39) 

0.00000 

(0.0000) 

0.01775 

(0.0628) 

0.07663 

(0.0805) 

0.18783 

(0.1017) 

0.36774 

(0.1186) 

0.64159 

(0.1448) 

1.05057 

(0.1771) 

1.66394 

(0.3435) 

2.54197 

(2.1235) 


Table 2f. Comparison of deflections calculated from SPAR with those calculated from different 


deflection equations; tapered cantilever tubular beam; c n Icq = 0.25/4; n = 8. 


Theories 

Jo 

(error, 

percent) 

Jl 

(error, 

percent) 

J2 

(error, 

percent) 

J3 

(error, 

percent) 

J4 

(error, 

percent) 

J5 

(error, 

percent) 

J6 

(error, 

percent) 

J7 

(error, 

percent) 

Jg 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01952 

(0.0000) 

0.08017 

(0.0000) 

0.19668 

(0.0000) 

0.38984 

(0.0000) 

0.69566 

(0.0000) 

1.18157 

(0.0000) 

2.00286 

(0.0000) 

3.58289 

(0.0000) 

Nonuniform, 
eq. (36) 

0.00000 

(0.0000) 

0.01793 

(0.0444) 

0.07826 

(0.0533) 

0.19445 

(0.0622) 

0.38766 

(0.0608) 

0.69340 

(0.0631) 

1.17872 

(0.1103) 

2.00100 

(0.0519) 

3.43492 

(4.1299) 

Slightly 
nonuniform, 
eq. (37) 

0.00000 

(0.0000) 

0.01785 

(0.0466) 

0.07771 

(0.0687) 

0.19276 

(0.1094) 

0.38359 

(0.1744) 

0.68447 

(0.3123) 

1.15879 

(0.6358) 

1.95019 

(1.4700) 

3.32381 

(7.2310) 

First-order, 
eq. (38) 

0.00000 

(0.0000) 

0.01789 

(0.0455) 

0.07796 

(0.0617) 

0.19354 

(0.1216) 

0.38546 

(0.1222) 

0.68850 

(0.1998) 

1.16756 

(0.3910) 

1.97136 

(0.8792) 

3.34045 

(6.7666) 

Second- 
order, 
eq. (39) 

0.00000 

(0.0000) 

0.01793 

(0.0444) 

0.07824 

(0.0539) 

0.19436 

(0.0648) 

0.38743 

(0.0673) 

0.69284 

(0.0787) 

1.17719 

(0.1222) 

1.99555 

(0.2040) 

3.40914 

(4.8494) 


Note from tables 2a-f that the deflections predicted from the four deflection equations (36)-(39) are 
fairly close. Note also that the SPAR program provided slightly larger deflections than the corresponding 
deflections predicted from the four deflection equations. The reason for this difference is that the SPAR 
deflections contain additional deflection components of transverse shear, which are neglected in the 
formulation of deflection equations (36)-(39) (see graphical demonstrations in reference 2, figures 17 
and 19). Also note that the prediction accuracy of deflection equations (36) and (39) are slightly better 
than the comparable accuracy of deflection equations (37) and (38). 

The deflection data of tables 2a-f calculated from deflection equations (36)-(40) are plotted in 
figures 5a-d, respectively, for different tapered cantilever tubular beams, and are compared with the 
corresponding deflection curves calculated from the SPAR program. The small prediction errors shown in 
tables 2a-f are almost inconspicuous for depth ratio range, 4 / 4 > c n / cq > 1 / 4 . For the cases of low 
depth ratios, c n i cq = (0.5/4, 0.25/4), however, the predicted deflection curves near the beam tip regions 
show slight deviations from the corresponding SPAR deflection curves because of the insufficient strain 
data points used. As discussed in the section entitled, “Domain Density,” the shape prediction accuracy 
near the beam tip regions could be greatly improved by increasing the domain density (increasing the 
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number of strain stations). Overall, the good graphical correlations between the predicted and SPAR- 
generated deflection curves provide confidence in the high accuracy of deflection equations (36)— (39) 
formulated for the nonuniform cantilever beam. 

Error Curves 

Figure 6 shows the beam tip prediction errors (tables 2a-f) based on different deflection equations 
plotted as functions of depth ratio, c n / cq . For the uniform beam case (c n /cg = 4/4), the deflections 
calculated from deflection equations (37)-(39) are identical, indicating the mathematical accuracy of 
deflection equations (37) — (39) at the limit of the uniform beam case. For the tapered beam cases, all the 
deflection equations (36)— (39) predicted reasonably close deflections (a prediction error of approximately 
1 percent) in the region of the depth ratios, 3/4 >c n /cg >2/4. As the beam depth ratio decreases from 
c n t c o =2/4, the prediction errors continue to increase and reach a respective maximum at the lowest 
depth ratio, c h /cq= 0.25/4. The beam tip shape prediction errors at c n / cq = 0.25/4 of both the 
nonuniform deflection equation (36) and second-order deflection equation (39) are comparable, and are 
nearly 60 percent of the prediction errors of both the slightly nonuniform deflection equation (37) and the 
first-order deflection equation (38). 


DOMAIN DENSITY 

As shown in figure 4, as the beam depth ratio, c n / cq , decreases, the strain curve continues to bend 
more and reaches a sharp bend near the beam tip region for the lowest depth ratio, c n /cg= 0.25/4. To 
obtain more accurate strain curves with a piecewise linear approximation (eq. (5)), the strain data points 
must be increased. Increasing the domain density, n, implies increasing the number (1 + n) of strain 
stations. Therefore, the domain density, n, was increased from n = 8 to n = 16 to obtain additional strain 
data points for generating more accurate strain curves, thereby achieving better shape predictions. 

Improved Strain Curves 

Figure 7 compares the SPAR-generated strain curves based on the domain densities of n = 8 and 
n = 16. For the depth ratios, c n /cg= 4/4-1/4, the effect of increasing the domain density, n, is practically 
inconspicuous. Increasing the value of n, however, has a marked effect on improving the strain curves 
near the beam tip regions especially for the lower depth cases, c n / cq = 0.5/4 and c n / cq = 0.25/4. Again, 
the nonzero strains at the free beam tip are caused by discretization of the beam (theoretically, the free 
beam tip strain should be zero). 


Improved Deflection Data 

As shown in figure 7, the beam deflections were recalculated from only deflection equations (36) 
and (39) (fig. 6) for different tapered tubular beams using the strains based on n = 16. Tables 3 a— f list the 
results, with prediction errors indicated in parentheses. 
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Table 3a. Comparison of deflections calculated from SPAR with those calculated from different 


deflection equations; uniform cantilever tubular beam; c n l cq = 4/4; n = 16. 


Theories 

Jo 

(error, 

percent) 

J2 

(error, 

percent) 

J4 

(error, 

percent) 

J6 

(error, 

percent) 

Js 

(error, 

percent) 

JlO 

(error, 

percent) 

y 12 

(error, 

percent) 

y 14 

(error, 

percent) 

y i6 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01787 

(0.0000) 

0.06340 

(0.0000) 

0.13368 

(0.0000) 

0.22450 

(0.0000) 

0.33181 

(0.0000) 

0.45174 

(0.0000) 

0.57933 

(0.0000) 

0.71162 

(0.0000) 

Nonuniform, 
eq. (36) 










Second- 
order, 
eq. (39) 

0.00000 

(0.0000) 

0.01576 

(0.2965) 

0.06034 

(0.4300) 

0.12962 

(0.5705) 

0.21948 

(0.7055) 

0.32580 

(0.8446) 

0.44448 

(1.0202) 

0.57139 

(1.1158) 

0.70243 

(1.2914) 


Table 3b. Comparison of deflections calculated from SPAR with those calculated from different 


deflection equations; tapered cantilever tubular beam; c n Icq = 3/4; n = 16. 


Theories 

Jo 

(error, 

percent) 

J2 

(error, 

percent) 

J4 

(error, 

percent) 

J6 

(error, 

percent) 

Js 

(error, 

percent) 

JlO 

(error, 

percent) 

J 12 
(error, 
percent) 

J 14 
(error, 
percent) 

J 16 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01823 

(0.0000) 

0.06697 

(0.0000) 

0.14572 

(0.0000) 

0.25204 

(0.0000) 

0.38294 

(0.0000) 

0.53428 

(0.0000) 

0.70092 

(0.0000) 

0.87634 

(0.0000) 

Nonuniform, 
eq. (36) 

0.00000 

(0.0000) 

0.01627 

(0.2237) 

0.06417 

(0.3195) 

0.14199 

(0.4256) 

0.24738 

(0.5318) 

0.37726 

(0.6482) 

0.52758 

(0.7645) 

0.69312 

(0.8901) 

0.86711 

(1.0532) 

Second- 
order, 
eq. (39) 

0.00000 

(0.0000) 

0.01627 

(0.2237) 

0.06417 

(0.3195) 

0.14199 

(0.4256) 

0.24738 

(0.5318) 

0.37726 

(0.6482) 

0.52759 

(0.7634) 

0.69313 

(0.8889) 

0.86712 

(1.0521) 


Table 3c. Comparison of deflections calculated from SPAR with those calculated from different 


deflection equations; tapered cantilever tubular beam; c n Icq = 2/4; n = 16. 


Theories 

Jo 

(error, 

percent) 

J2 

(error, 

percent) 

J4 

(error, 

percent) 

J6 

(error, 

percent) 

Js 

(error, 

percent) 

JlO 

(error, 

percent) 

J 12 
(error, 
percent) 

J 14 
(error, 
percent) 

J 16 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01864 

(0.0000) 

0.07113 

(0.0000) 

0.16057 

(0.0000) 

0.28831 

(0.0000) 

0.45510 

(0.0000) 

0.65965 

(0.0000) 

0.89759 

(0.0000) 

1.15801 

(0.0000) 

Nonuniform, 
eq. (36) 

0.00000 

(0.0000) 

0.01682 

(0.1572) 

0.06859 

(0.2193) 

0.15719 

(0.2919) 

0.28407 

(0.3661) 

0.44983 

(0.4551) 

0.65325 

(0.5527) 

0.88981 

(0.6718) 

1.14828 

(0.8402) 

Second- 
order, 
eq. (39) 

0.00000 

(0.0000) 

0.01682 

(0.1572) 

0.06859 

(0.2193) 

0.15719 

(0.2919) 

0.28407 

(0.3661) 

0.44983 

(0.4551) 

0.65324 

(0.5535) 

0.88980 

(0.6727) 

1.14826 

(0.8420) 


Table 3d. Comparison of deflections calculated from SPAR with those calculated from different 


deflection equations; tapered cantilever tubular beam; c n Icq = 1/4; n = 16. 


Theories 

Jo 

(error, 

percent) 

J2 

(error, 

percent) 

J4 

(error, 

percent) 

J6 

(error, 

percent) 

Js 

(error, 

percent) 

JlO 

(error, 

percent) 

J 12 
(error, 
percent) 

J 14 
(error, 
percent) 

J 16 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01912 

(0.0000) 

0.07598 

(0.0000) 

0.17923 

(0.0000) 

0.33816 

(0.0000) 

0.56506 

(0.0000) 

0.87454 

(0.0000) 

1.28196 

(0.0000) 

1.78417 

(0.0000) 

Nonuniform, 
eq. (36) 

0.00000 

(0.0000) 

0.01739 

(0.0970) 

0.07364 

(0.1312) 

0.17612 

(0.1743) 

0.33427 

(0.2180) 

0.56011 

(0.2774) 

0.86835 

(0.3469) 

1.27375 

(0.4602) 

1.77114 

(0.7303) 

Second- 
order, 
eq. (39) 

0.00000 

(0.0000) 

0.01739 

(0.0970) 

0.07364 

(0.1312) 

0.17611 

(0.1749) 

0.33426 

(0.2186) 

0.56009 

(0.2786) 

0.86829 

(0.3503) 

1.27365 

(0.4658) 

1.77124 

(0.7247) 
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Table 3e. Comparison of deflections calculated from SPAR with those calculated from different 


deflection equations; tapered cantilever tubular beam; c n Icq = 0.5/4; n = 16. 


Theories 

Jo 

(error, 

percent) 

J2 

(error, 

percent) 

J4 

(error, 

percent) 

J6 

(error, 

percent) 

Js 

(error, 

percent) 

JlO 

(error, 

percent) 

y 12 

(error, 

percent) 

y 14 

(error, 

percent) 

y i6 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01938 

(0.0000) 

0.07872 

(0.0000) 

0.19047 

(0.0000) 

0.37082 

(0.0000) 

0.64535 

(0.0000) 

1.05517 

(0.0000) 

1.67292 

(0.0000) 

2.59712 

(0.0000) 

Nonuniform, 
eq. (36) 

0.00000 

(0.0000) 

0.01770 

(0.0647) 

0.07649 

(0.0859) 

0.18752 

(0.1136) 

0.36721 

(0.1390) 

0.64077 

(0.1763) 

1.04963 

(0.2133) 

1.66517 

(0.2984) 

2.57675 

(0.7843) 

Second- 
order, 
eq. (39) 

0.00000 

(0.0000) 

0.01770 

(0.0647) 

0.07648 

(0.0862) 

0.18751 

(0.1140) 

0.36719 

(0.1398) 

0.64071 

(0.1787) 

1.04949 

(0.2187) 

1.66478 

(0.3134) 

2.57544 

(0.8348) 


Table 3f. Comparison of deflections calculated from SPAR with those calculated from different 


deflection equations; tapered cantilever tubular beam; c n Icq = 0.25/4; n = 16. 


Theories 

Jo 

(error, 

percent) 

J2 

(error, 

percent) 

J4 

(error, 

percent) 

J6 

(error, 

percent) 

Js 

(error, 

percent) 

JlO 

(error, 

percent) 

J 12 
(error, 
percent) 

J 14 
(error, 
percent) 

J 16 

(error, 

percent) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01952 

(0.0000) 

0.08017 

(0.0000) 

0.19668 

(0.0000) 

0.38984 

(0.0000) 

0.69566 

(0.0000) 

1.18157 

(0.0000) 

2.00286 

(0.0000) 

3.58289 

(0.0000) 

Nonuniform, 
eq. (36) 

0.00000 

(0.0000) 

0.01786 

(0.0463) 

0.07802 

(0.0600) 

0.19387 

(0.0784) 

0.38651 

(0.0929) 

0.69143 

(0.1181) 

1.17563 

(0.1658) 

1.99425 

(0.2403) 

3.53790 

(1.2557) 

Second- 
order, 
eq. (39) 

0.00000 

(0.0000) 

0.01786 

(0.0463) 

0.07801 

(0.0603) 

0.19385 

(0.0790) 

0.38648 

(0.0938) 

0.69135 

(0.1203) 

1.17539 

(0.1725) 

1.99331 

(0.2665) 

3.53116 

(1.4438) 


The deflection data shown in tables 3a-f are plotted in figure 8. Notice from figures 8a and 8b that the 
beam tip prediction errors using n = 16 did indeed decrease from those using n = 8 (tables 2a-f)), 
especially for the lower depth ratio cases, c n / cq= 0.25/4, 0.5/4. 


Error Comparisons 

Tables 4a-f compare the beam tip prediction errors based on n = 8 and n = 16. 


Table 4a. Comparison of beam tip prediction errors based on different deflection equations using strain 
station densities (n = 8, 16); tapered cantilever tubular beam; c n Icq = 4/4. 




n = 8 


n = 1 6 


Theories 

y n > in - 

Error, percent 

y n > in - 

Error, percent 

Error, percent of 
n = 8 

SPAR 

0.71162 

0.00 

0.71162 

0.00 

0 

Nonuniform, 
eq. (36) 

— 

— 

— 

— 

— 

Second-order, 
eq. (39) 

0.70230 

1.31 

0.70243 

1.29 

98 
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Table 4b. Comparison of beam tip prediction errors based on different deflection equations using strain 


station densities ( n = 8, 16); tapered cantilever tubular beam; c n Icq 

= 3/4. 



n = 8 


n =16 


Theories y 5 i n . 

Error, percent 

y n ’ in ' 

Error, percent 

Error, percent of 
n = 8 

SPAR 0.87634 

0.00 

0.87634 

0.00 

0 

Nonuniform, 0.86568 

1.22 

0.86711 

1.05 

86 

eq. (36) 





Second-order, 0.86567 

1.22 

0.86712 

1.05 

86 

eq. (39) 





Table 4c. Comparison of beam tip prediction errors based on different deflection equations using strain 

station densities ( n = 8, 16); tapered cantilever tubular beam; c n Icq 

= 2/4. 



n - 8 


n = 16 


Theories y^ 5 i n . 

Error, 

y n , in- 

Error, percent 

Error, percent of 


percent 



n = 8 

SPAR 1.15801 

0.00 

1.15801 

0.00 

0 

Nonuniform, 1.14592 

1.04 

1.14828 

0.84 

81 

eq. (36) 





Second-order, 1.14518 

1.1 

1.14827 

0.84 

76 

eq. (39) 





Table 4d. Comparison of beam tip prediction errors based on different deflection equations using strain 

station densities ( n = 8, 16); tapered cantilever tubular beam; c n Icq 

= 1/4. 



n = 8 


n = 16 


Theories y^ 5 i n . 

Error, 

y n , in. 

Error, percent 

Error, percent of 


percent 



n = 8 

SPAR 1.78417 

0.00 

1.78417 

0.00 

0 

Nonuniform, 1.76228 

1.23 

1.77144 

0.71 

58 

eq. (36) 





Second-order, 1.76105 

1.30 

1.77124 

0.72 

55 

eq. (39) 






Table 4e. Comparison of beam tip prediction errors based on different deflection equations using strain 
station densities ( n = 8, 16); tapered cantilever tubular beam; c n Icq = 0.5/4. 




n = 8 


n = 16 


Theories 

J„,in. 

Error, 

percent 

y n > in - 

Error, percent 

Error, percent of 
n = 8 

SPAR 

2.59712 

0.00 

2.59712 

0.00 

0 

Nonuniform, 
eq. (36) 

2.54871 

1.86 

2.57675 

0.78 

42 

Second-order, 
eq. (39) 

2.54197 

2.12 

2.57544 

0.83 

39 


23 











Table 4f. Comparison of beam tip prediction errors based on different deflection equations using strain 
station densities (« = 8, 16); tapered cantilever tubular beam; c n Icq = 0.25/4. 




n = 8 


n = 16 


Theories 

y n > in - 

Error, 

percent 

T„,in. 

Error, percent 

Error, percent of 
n = 8 

SPAR 

3.58289 

0.00 

3.58289 

0.00 

0 

Nonuniform, 
eq. (36) 

3.43492 

4.13 

3.53790 

1.26 

31 

Second-order, 
eq. (39) 

3.40914 

4.85 

3.52116 

1.44 

30 


Figure 9 plots the beam tip deflection prediction errors based on n = 8 and n = 16 (tables 4a-f) for 
comparison. Note that the effect of increasing the domain density, n, on the reduction of prediction errors 
becomes more conspicuous as the depth ratio, c n / cq , decreases. At the lowest depth ratio, 
c n Icq= 0.25/4, the beam tip prediction errors based on n = 16 were reduced to almost 30 percent of those 
based on n = 8. 


DISCUSSION 

For the present shape predictions of the tapered cantilever tubular beam, all four displacement transfer 
functions (deflection equations (36)-(39)) were found to provide quite accurate shape predictions. The 
slightly nonuniform deflection equation (37) and the first-order deflection equation (38) have comparable 
prediction accuracy and are recommended for shape predictions of slightly nonuniform beam cases. For 
shape predictions of nonuniform beam structures, both the exact deflection equation (36) and second- 
order deflection equation (39) can be used, because both deflection equations can provide comparable 
shape prediction accuracy. 

Note that when the fiber-optic strain-sensing system is used, off-the-shelf sensor density is high, 
thereby improving shape prediction accuracy. This high sensor density can be decimated by an integer 
factor to reduce computational requirements while sacrificing prediction accuracy. Thus, for improving 
shape prediction accuracy, the previously specified number {n + 1) of strain data extraction points can be 
increased as desired, up to a level corresponding to the off-the-shelf sensor density, by means of a simple 
command. By increasing the domain density, n, the assumed piecewise-linear strain curves will approach 
the associated true smooth strain curves (especially for the highly tapered beam tip region where the 
bending strains reach a maximum and then drop rapidly toward the beam tip), and thus greatly improve 
the shape prediction accuracy. 


SUMMARY 

The first- and the second-order displacement transfer functions were formulated for the nonuniform 
cantilever beam through binomial series expansions of the beam depth function. The principal results are 
as follows: 


1 . Through binomial series expansions of the axially varying beam depth function, the mathematical 
indeterminacy problem at the limit case of a uniform beam was successfully eliminated. 

2. The deflections calculated from all four displacement transfer functions were reasonably close and 
were slightly lower than the corresponding deflections generated by the Structural Performance 
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And Resizing (SPAR) program. The SPAR-generated deflections contained additional 
components induced by the transverse shear effect, which was ignored in the analysis. 

3. For the uniform beam case (c w /cq = 4/4), the deflections predicted from different displacement 
transfer functions were identical, confirming the mathematical accuracy of all the displacement 
transfer functions formulated. 

4. The shape prediction accuracy of the nonuniform and second-order displacement transfer functions 
were slightly better than the comparable accuracy of both the slightly nonuniform and first-order 
displacement transfer functions. 

5. For the tapered beam shape predictions, using the second-order displacement transfer function 
instead of the first-order deflection equation did not change the beam tip prediction error for the 
uniform beam (zero taper angle). For the tapered beam, however, the prediction error continued to 
decrease with increasing taper angle, reaching a 28-percent error reduction at the peak taper angle. 

6. For the tapered beam shape predictions, using the nonuniform displacement transfer function 
instead of the slightly nonuniform deflection equation decreased the beam tip prediction errors 
with increasing taper angle, reaching a 43-percent error reduction at the peak taper angle. 

7. For the tapered beam, the shape prediction accuracy could be improved by increasing the domain 
density. By doubling the domain density, one could reduce the average beam tip prediction errors 
by 2 percent at a zero taper angle (uniform beam). The error reduction continued to improve with 
increasing taper angle, reaching a 69-percent error reduction at the peak taper angle. 
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FIGURES 



Figure 1. Nonuniform cantilever beam instrumented with equally spaced surface bending strain sensors. 


P= 100 lb' 



/ =100.5 in. Nodes 3673 

t = 0.02296 in. Four-node elements 3636 

c 0 = 4 in. Three-node elements 36 

c n = 4, 3, 2, 1, 0.5, 0.25 in. 110283 

Figure 2. Finite-element model for tapered cantilever tubular beam subjected to a tip vertical load of 
P= 1001b. 
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(a) c n/ c 0 


= 4/4 


110284 



(b) c n /c 0 = 3/4 


110285 



(C) c n /c 0 = 2/4 


110286 



Figure 3. Deformed shapes, generated from the Structural Performance And Resizing (SPAR) program, 
of different tapered cantilever tubular beams; P = 100 lb. 
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P= 100 lb 



(f) c„/c 0 = 0.25/4 
Figure 3. Concluded. 



Figure 4. Surface bending strains calculated from the Structural Performance And Resizing (SPAR) 
element stress outputs for different tapered cantilever tubular beams; n = 8; P= 100 lb. 
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(a) Deflections calculated from deflection equation (36), for nonuniform beams, compared with those 
calculated from SPAR (deflection equation (40) was used for the uniform beam). 



(b) Deflections calculated from deflection equation (37), for slightly nonuniform beams, compared with 
those calculated from SPAR (deflection equation (40) was used for the uniform beam). 

Figure 5. Comparison of deflections calculated from deflection equations with those calculated from the 
Structural Performance And Resizing (SPAR) program for different tapered cantilever tubular beams; 
« = 8; R = 100 lb. 
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(c) Deflections calculated from first-order deflection equation (38) compared with those calculated from 
SPAR (deflection equation (40) was used for the uniform beam). 



(d) Deflections calculated from second-order deflection equation (39) compared with those calculated 
from SPAR (deflection equation (40) was used for the uniform beam). 

Figure 5. Concluded. 
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Figure 6. Beam tip deflection prediction errors associated with different displacement transfer functions 
plotted as functions of depth ratio, c n / Cq . 


Strain station, ; 



Strain station, i 


110296 


Figure 7. Comparison of bending strain curves for different tapered cantilever tubular beams based on 
n = 8 and n = 16; P = 100 lb. 
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Strain station, i 

110297 

(a) Deflections calculated from deflection equation (36), for nonuniform beams, compared with those 
calculated from SPAR (deflection equation (40) was used for the uniform beam). 



Strain station, i 

110298 

(b) Deflections calculated from second-order deflection equation (39) compared with those calculated 
from SPAR. 

Figure 8. Comparison of deflections calculated from deflection equations with those calculated from the 
Structural Performance And Resizing (SPAR) program for different tapered cantilever tubular beams; 
n= 16; P= 1001b. 


32 



Depth ratio, cjc 0 

' 110299 


Figure 9. Comparison of beam tip deflection prediction error curves associated with different 
displacement functions based on n = 8 and n = 16. 
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APPENDIX A 


DERIVATIONS OF THE FIRST-ORDER SLOPE 
AND DEFLECTION EQUATIONS 


Appendix A presents the detailed mathematical derivations of the first-order slope equations and 
first-order deflection equations for the nonuniform cantilever beam. 


Basic Slope Equations 

The slope, tan0(x)(= dy / dx) , of the nonuniform beam at the axial location, x, within the domain, 
x ; _| <x<x- (fig. 1), can be obtained by integrating the beam differential equation (3) once, and 
enforcing the continuity of the slope at the inboard strain station, jc -_j , as 


tan 0(x) = j 

J X: 


,2 

d y 


dx 


, x i~ 1 dx A 
Integration ofeq. (3) 


+ tan0^_i = J* 

Slope at x,-_ j 


£(*) 
x i - 1 c(x) 


Slope increment 
above tan Q t _\ 


dx + tan0 ; _| ; (x,_j < x < Xj) 

Slope at x,-_i 


(Al) 


in which tan 0._j is the slope at the inboard strain station, x-_j . 

In light of the following linear assumptions of both c(x) and £(x) within the domain, x ; _| < x < x. , 


Ci) X ~ Xi - 1 

Al 

; (x ; -_i <x<Xj) 

(A2) 

, x - Xl_ | 

£;) — L 

A/ 

; (x-_, <x<x.) 

(A3) 


equation (Al) becomes 


1- 


f £■ A 
1-^ 


tan 


(Kx) = *±\ x 

r- i •'I; 


V £ i - 1 


X-Xj-l 
Al 


Ci - 1 J A-1 


1- 


f c- ^ 
v c i-i ; 


Al 


- r/x + tan 0,_] ; (x •_ , < x < x • ) 


(A4) 
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First-Order Slope Equations 

Because the coefficient, 1 — (c- / c -_^ ) , of the slope term in the denominator of equation (A4) is 


considered small, one can expand the factor, 
first- order term as 


i-l 


1- 


1--^- 


x-x 


i-l 


V L i-1 J 


M 


in binomial series up to the 


tan0(x) = [ 

r i - 1 ix i - 1 


1- 


1--^- 




x — x 


i-l 


_ e,-_ i f.v 


f 

i-l Jx ‘- 1 


1 + 


1--^- 


v fc /-i ; 

A f 


v c i- 1 ; 


1--^- 


A/ 

Y 


1 + 


1--^- 


V C i-U 


x — x. 


i-l 


K £ i- 1 7 


x-x 


i-l 


A/ 


1--^- 


A/ 

V 


+ ... 


t/x + tan0 


V £ /-lA l i-l) 




C; 


( X ~ X i-l) 


i-l 


+ .... 


(A/)' 


</x + tan0._| (A5) 


After the integration is carried out, equation (A5) becomes 


tan 0(x) = 


£ i-l 


c; 


i—l 


(x-x,_ i) + 


c; 


v fy — 1 c i — 1 / 


(X-X;_i) 


2A/ 


1 -- 


1- 


c; 


(X-X;_i) J 


V c i-l7 


3(A iy 


+ tan 0/_ 


i—l 


(x._, <X<X.) (A6) 

which is the first-order slope equation for the domain, x-_j < x < x- . 

At the strain station, x = x ; - , we have (x - -x-_j) = A l , and equation (A6) becomes 

A/e f _ i r 


tan 9i = tan 0(x,- ) = 


c i-l 
A / 


i + i 

2 




£ i c i 




6c, ■ _ 


i—l 


V £ i — 1 c i — 1 J 

( 


1 £± - 


V £ i-l 7 v c i-l 7 


1 — Cj - 


+ tan 0,_ 


i—l 


6-3— 2 


c i-l 


1 — 

V c i-l 7 


£ i - 1 + 


3 + 2 


-iA 

J 


c 

Y 

1 1 

1-- 

£ i 

£ . 1 

l 1 

c i- lj. 

J J 


A-l 


(i = 1, 2, 3, ..., m) (A7) 


After the terms are grouped, equation (A7) becomes 
A/ 


tan 6- = 

1 6c 


i-l 


( c i ^ 
4 E 

l C i-D 


£ i - 1 + 


' c. ' 
5-2 — - 


V 


'i-l 7 


+ tan0 ; _ 1 ; (i = l, 2, 3, ..., «) 


(A8) 
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Equation (A8) is equation (17) in the main section of this report, and the first-order slope equation is 
written in descending recursion form. Applying the descending recursion relationships, one can write 
equation (A8) in summation form as 


„ A/ J, 1 
tanO, = — £ 


6 Pl C J~ 1 


4- 


c ■ 


7-1 


£ 7-i + 


5-2 


C j 


-j - 1 


tan 0 q 

=0 for canti- 
lever beams 


; ()'= 1, 2, 3, ..., n) 


(A9) 


Equation (A9) is equation (18) in the main section of this report, and the first-order slope equation is 
written in series summation form. 

For the uniform beam ( c ; - = c-_ j = c ), equation (A9) degenerates into the form 


A / 


tan 0j = — ^ (£ /_i + £ / ) + tan0 o ; (z = 1, 2, 3, ..., n) 


2c 


7=1 


=0 for canti- 
lever beams 


(A 10) 


which is identical to equation (25) of reference 2. 

First-Order Deflection Equations 

The deflection, y(x) , at the axial location, x, within the domain, x-_^ < x < x ; - (fig. 1), can be 
obtained by integrating equation (A6) and enforcing the continuity of the deflection at the inboard strain 
station, x-_j . Namely, 


rX 

y(x ) = I tan 6{x)dx + y-_ | 

Jx i - r v ' 


Eq. (A6) 


_ g /-l f* 

. J ' 


i - 1 »-l 


(x-x-_|) + 


£ i c / 


V £ /-l c i-l 7 


( x ~ x i- 1> 


2 A/ 


1--^- 


V £ i-l 7 


1--^- 




^ W-i 7 


dx 


3(A/) 

+ J* tan 0 ; _, 7x + y._, 

(x/—| <x<Xj) (All) 


h-i 
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in which y ( — ] is the deflection at the inboard strain station, jc-_j . After the integration is carried out, 
equation (All) becomes 


y(x) = 


£ i - 1 


12c 


i - 1 


6(X~ + 2 


£ i 


V £ /- 1 c i-\) 


( x ~ x i- l> 


A/ 


1 — £j - 


~"i- 1 y 


1--^- 


~i- 1 y 


( x - x f-|) 


(A/)" 


+ (x-x._ 1 )tan0._ 1 + j._ 1 

0 ; -_l <x<x j ) (A 12) 


Equation (A 12) is the first-order deflection equation for the domain, x-_j <x<Xj. 

At the strain station, x = x- , we have (x- — x-_j) = A/ , and equation (A 12) takes on the form 


yi = y( x i)= 


(A /re M 


12c ; ._i 


£; C- 

6 + 2 — 2 — 

£ /'-l c /-l 


2 £ / c i | c i £ i 

v £ /-i c /-i c /-i £ /-i y 


+ _y -_| + A/ tan0 ; ._ 


i— 1 


(f = l, 2, 3, ..., «) (A13) 


After the terms are grouped, equation (A13) becomes 


yi = 


(A/)" 


1 


c; 


5-- 
V c i-i ) 


( 


£ i - 1 + 


3-AA 

v c i-i ; 


+ y,_l + A/tan0 ; -_i ; (/' = 1, 2, 3, ..., «) (A 14) 


Equation (A 14) is equation (21) in the main section of this report, and the first-order deflection equation 
is written in descending recursion form. 


For the uniform beam ( c ; - = c •_ j = c ), equation (A14) degenerates into 


>7 = < AT ) -(2 £ M + £ /) + >7-1 + A/ tan ; (i = 1, 2, 3, ..., n) (A15) 

6c 

which is identical to equation (26) of reference 2. 

Final Form of Equation (A14) 

Substituting the slope equation (A8) into deflection equation (A 14), and applying the descending 
indicial recursion relationships, one can write the resulting deflection equation for each index, 
i(= 1, 2, 3, ..., n), as follows: 
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For i=l, 


(M) z 1 cj 


31 12 c 0 l c o J ° 


— 5 £q + 3 £ 1 +>’o + ^ tan ®o 


For i = 2, 


(A/) 2 if c 2 ^ f c 2 ^ (A/) 2 if c, ) f c, ) 

y 2 -——— 5 ~— £, + 3 -— £ 2 + — — 5 “— £ o + 3_ ~ £ l +v O + A/tan0 O 

u c i l c \ ) V c i) 1Z c o { c o ) \ c o J 


6 c 0 l c 0 


(Mr l c, c, 

H 4 £ n + 5-2 — £, + A / tan 0, 


A/ tan 0 , 


After the terms are grouped, equation (A 17) becomes 


(A/r i L c 2 


12 ^ Cj 


c lJ c o l c 0 


C C \ c cl 

5 — — £j + 3 — — £ 0 -4 13-3— £ 0 + 13-5— £j t +v 0 + 2A/tan0 o (A18) 


For i = 3, 


(A/r i c 3 


>3 12 c 2 l c 2 2 


— 5 £o + 3 £ 


, . „2 , iy a f t i (1 „2 , r/ A ( \ 

(Air 1 c 7 c 7 (Al) 1 c i c i 

+ — 5 — ± £,+ 3 — — £ 7 + — 13-3— £ n + 13-5— £, + v n + 2 A/ tan 0, 

12 C, Cl 1 c, 2 12 Cn Cn 0 C n 1 ' 0 * 


1 + >’0 


(a iy if. c 2 


6 C 1 l C 1 


(A/)" 1 ci 


+ — 4 — — £, + 5-2— £o + — 4 £n+ 5-2— £i +A/tan0, 


6 c oLl c 0 

Al tan 0 ^ 


After the terms are grouped, equation (A 19) becomes 


< A ') Z 4f 5 _fi) e , + (3-fil 3 l + A^i[f,3-3^) ei + fl3-^] e 

Co Co Co J 12 c, c, c. 


>3 12 c 2 l c 2 J 2 


. ..2 1 IY A f A 

+ — 21-5— £ 0 + 23-9— £j +v o + 3A/tan0 o 

1 2 Cn Cn Cn 
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For i = 4, 


>’4 


(A l) 2 1 


+ 


12 c 3 
(A/) 2 1 


17 5 


c 

\ 

“ 

5— — 

£ 3 + 

3- 

O 

e 4 

A C 3 J 


V 

c 3y 



12 


17 5 

Cn 

5 — A 

£>■) + 

( 

3- 

0 

£9 

, (A/) 2 1 " 

y x 

Cn 

13-3 — 

£ 1 + 

7 x 

Cn 

13-5 — 

£9 

A C 2 J 

Z 

V 

c 2y 

J 

12 c 3 

l c lJ 

l c lJ 

z 


(A/) 2 1 


12 Cn 


\( \ 


7 

\ 


Cl 

21-5 — 

£ 0 + 

23- 

c, 

- 9 — 

£ 1 

A c o J 

1 

c 0 J 



+ Vq + 3 A / tan Oq 


^3 


+ 


(A/) 2 1 


6 Co 


17 \ 


( \ 

- 
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y 
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( \ 

~ 

Cn 

4 — A 

e 2 + 

Cn 

5-2 — 

£ 3 
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4- 

O 

£ 1 + 

Cn 

5-2 — 

£ 2 

A C 2 J 


l c 2 J 


6 c 3 

A 

C l7 


l c lJ 



1 

v 

4 — 

A 

0. 

c* _i_ 

( \ 
Cl 

S— ? — L 


c 0 . 

‘-r 

A 

c 0 j 

£ 0 + 

J z, 

l c 0 J 

£ 1 


+ A/ tan 8, 


0 


A / tan On 


After the terms are grouped, equation (A21) becomes 


(A21) 


>’4 = 


(A/) 2 1 


12 c 3 
(A/) 2 1 


12 Cn 


5 — — 
c - 


13-3 — 


\ 

7 


“ 

£ 3 + 

3- 

O 

e 4 

7 

V 

c 3y 



L 2 7 


(A/) 2 1 


12 c, 


e 2 + 
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( c A 

13-5— 

V C 2) 

. 5 


(A/) 2 1 


29-7 — 

c 0j 


V 


7 


£ 0 + 


12 c, 


A 


21-5 — 
c, 


7 


£ 1 + 


\ 




33-13 — 

c 0 ) 


V 


17 

+ Jq +4 A/ tan0Q 


23-9 — 

C 1 7 


(A22) 


Observing the indicial behavior in equations (A 16), (A18), (A20), and (A22), one can write the 
generalized first-order deflection equation for any index, z(= 1, 2, 3,..., n ) , in final summation form as 


y,-^i ' 


12 _ , c i-i 

7=1 i J 


5 + 8(y-l) -(l + 2(y-l)) 


- 7 - 7+1 


'i~J 


£ i~j 


3 + 100-1) -(1 + 40-1)) 


• 7 — 7+1 
C *— 7 




*—7+1 


+ Vq + (z)A/tan0Q 
=0 for cantilever beams 

0 = 1, 2, 3, ..., n) 


(A23) 


which is equation (23) in the main section of this report (called the first-order displacement transfer 
function). 
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For the uniform beam ( c- = c-_j = c ), equation (A23) degenerates into 


y ,- 


( a if 

12 c 


X {[ 5 + 80'-l)-(l + 2(y-l))]e f _ y . 

7=1 

+ [3 + 10(y - 1) - (l + 4(y - 1))] £ ; _ /+| } + j Q + (i)A/tan0 Q 


=0 for cantilever beams 

0 = 1, 2, 3, n) 

After the terms are grouped, the deflection equation (A24) for the uniform beam takes on the form 


(A 1)2 ^ 

yi = ~ 6 ^~ r; 

7=1 


(3y-l)e I ._ y .+(3y-2)e / _ l , 


1-7+1 


+ jVq +(/)A/tan0Q 
=0 for cantilever beams 

0 = 1, 2, 3, n) 


which agrees with equation (27) of reference 2. 

One can write out equation (A25) for typical indices, i (= 1, 2, 3, 4, ..., «), as follows: 
For i= 1, 


F] 


(A/)" 

6c 


(2£ 0 + £j ) + Vq + A/ tan 6q 


For i = 2, 


F 2 


= (A/y 

6c 

(A/f 

6c 


( 2£| + £-, + 5£q + 4£j ) + Vq + 2 A / tan 0 q 


■ (5£q + 6£| + £ 2 ) + Vq + 2 A l tan0Q 


For i = 3, 


F 3 


= (a iy 

6c 

(A l) 2 
6c 


- (2£q + £^ + 5£j + 4£t + 8£q + 7£j) + Vq + 3A l tan0Q 
[8£q + 6(2£j + £^ ) + £^ J + jvq + 3 A l tan0Q 


(A24) 


(A25) 


(A26) 


{All) 


(A28) 
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For i = 4, 


+ £4 + 5&2 + 4£^ + 8£^ + 7 £2 + 1 1£q + 1 0£^ ) + )?q + 4 A / tan 0q 

(A29) 

+ 6(3£j + 2&, + ) + £^ J + Vq + 4 A / tan 0 q 

Observing the functional behavior of equations (A26)-(A29) with an increasing index, i, one can write 
the generalized deflection equation for any index, /'(= 1, 2, 3, ..., n) , for the uniform cantilever beam as 

+ v 0 + O') A/ tan0 Q ; (i = 1, 2, 3, n) (A30) 

=0 for cantilever beams 

which agrees with equation (28) of reference 2. 


(3/ — 1)£q (' _ -/) £ y + £ / 

7=1 
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APPENDIX B 


DERIVATIONS OF THE SECOND-ORDER SLOPE 
AND DEFLECTION EQUATIONS 

Appendix B presents the detailed mathematical derivations of the second-order slope equations and 
second-order deflection equations for the nonuniform cantilever beam. 


Expanding the factor, 
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Integration of equation (Bl) results in 
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which is the general second-order slope equation for the domain, x/_] < x < x,- . 
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At the strain station, x = Xj , we have (xj - x/_] ) = A / , and equation (B2) becomes 
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After the terms are grouped, equation (B3) becomes 


tan 6; = 


M 


1 12c 


i-l 


4 „ 32 


9 — 4 — - — l- 
c i-l 


y c i-u 


£ i-i + 


4 „ 32 


13-10 — ^— + 3 
c i - 1 


V c i-l ) 


+ tan 0; 


i-l 


(i = 1, 2, 3, ..., ri) (B4) 


which is equation (28) in the main section of this report, the second-order slope equation in descending 
recursion form. Applying the descending recursion relationships, one can write equation (B4) in 
summation form as 
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which is the second-order slope equation in series summation form. 


43 



For the uniform beam, c,- = c;_i = c , equation (B5) degenerates into 


M yL 

tan Oj = — 2 j (£y_ i + £j)+ tan0 Q ; 0 = 1, 2, 3, n) 

=0 for canti- 
lever beams 


(B6) 


which is identical to equation (25) of reference 2. 

Second-Order Deflection Equations 

The second-order deflection equation is obtained by carrying out the integration of the second-order 
slope equation (B2) as 
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After integration of equation (B7), one obtains 
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which is the general second-order deflection equation for the domain, x/_ ] < x < x; . 
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At the strain station, x = X ( , we have (xj - x/_] ) = A / , and (B8) becomes 
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After the terms are grouped, equation (B9) becomes 
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which is equation (32) in the main section of this report, the second-order deflection equation in 
descending recursion form. 

For the uniform beam, q = = c , equation (B10) degenerates into 


y t = (2£ f _ | +£;) + A/ _| + A/ tan ; (i = 1, 2, 3, ..., w) ( B11 ) 

which agrees with equation (26) of reference 2. 

Final Form of Equation (B10) 

Substituting the slope equation (B5) into the deflection equation (BIO), and applying the descending 
indicial recursion relationships, one can write out the deflection equation (B10) explicitly for each index, 
/(= 1, 2, 3, ..., n), as follows: 


45 



For i= 1, 
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For i = 2, 
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After the terms are grouped, equation (B13) can be written as follows: 
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For i = 3, 


•^3 = 


(MY 

60 Co 


' 


7 3 

2 



7 3 

2 " 

' 


C-3 

27-9 — + 2 



e 2 + 

C-3 

18- 1 1 — + 3 



£ 3 


c 2 

v c 2 . 



c 2 

v c 2 . 




+ 


(MY 

60 c, 


27 - 9 — + 2 


-1 


7_ 32 


V C 1 7 


£ 1 + 


18- 1 1 — + 3 
C 1 


7„ 32 


V C 1 7 


+ 


(A/)- 

60c n 



j 

L 


J 

J 


f \2 1 



7 \ 

72-29 — + 7 



£ 0 + 

83-61 — + 18 

Yi 

c 0 

v c 0 y 


c 0 

v C 0 , 


+ 


Jq + 2 A/ tan0Q 


^2 


+ 


(MY 

12 c, 


' 


7 3 

2 ’ 



7 3 

2 ’ 

' 


9-4^- + 

^2 


£ 1 + 

c ? 

13- 10 — + 3 

f 2 


£ 2 


C 1 

v C l, 



C 1 

v C l, 




+ 


(A/)- 

12Cn 



J 

L 



J 

J 


7 3 

2 " 



7 3 

2 ’ 

c, 

9_4_L + 

f]_ 


£ o + 

13-10 — + 3 

Yl 


c 0 

l c o J 

_ 


1 

l c 0 J 

_ 


+ A/ tan 0, 


0 


Al tan 9~, 


After the terms are grouped, equation (B 15) can be written as follows: 
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(B 1 6 ) 
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For i = 4, 
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After the terms are grouped, equation (B17) becomes 
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Observing the indicial progressions in equations (B12), (B14), (B16), and (B18), one can write the 
generalized second-order deflection equation for any index, i{= 1, 2, 3, ..., n ) , in series summation form 
as 
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which is equation (34) in the main section of this report, called the second-order displacement transfer 
function. 
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For the uniform beam, c,- = c;_i = c, equation (B19) degenerates into 
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which can be written in a compact form as 
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which is identical to equation (A25) of Appendix A. 
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